Structure formation by cosmic strings with a cosmological constant 
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We investigate cosmic string models for structure formation which include a non-zero cosmological 
constant. We find that the background evolution of density perturbations and modifications to the 
scaling behaviour of the strings both act to increase the amount of power present on 100/i~^Mpc 
scales. We estimate the size of this effect using an analytic model for the evolution of the string 
network and find that a bias & ~ 2 can give acceptable agreement with the current observations. An 
interesting by product of these modifications is a broad peak in the cosmic microwave background 
(CMB) angular power spectrum around I = 400 — 600. 



PACS Numbers : 98.80. Cq, 95.35+d 

Topological defect theories, such as cosmic strings [Q, 
have for many years represented the only serious alterna- 
tives to the standard inflationary paradigm as theories for 
structure formation. In inflationary theories, structure 
in the universe is seeded by causal processes operating 
outside the framework of the Standard Big Bang (SBB) 
cosmology, which imprint a primordial spectrum of adi- 
abatic density perturbations on all scales. By contrast, 
in defect theories the universe is taken to be initially 
smooth and perturbations must be seeded by the causal 
action of an evolving network of defect seeds, within the 
background of the SSB cosmology. 

Recent work ||-^ has pointed to a serious problem for 
defect theories reconciling the amplitude of large scale 
CMB anisotropics with that in the matter distribution 
on scales around lOO/i^^Mpc, in a flat universe where 
the matter density is critical, that is il„i — 1. In refs. 

this problem was quantifled in terms of foioo, where 
bji is the bias between the cold dark matter (CDM) and 
the galaxy distribution in a sphere of radius i?/i^^Mpc 
that we must postulate in order to make a given theory 
consistent with the data. It was argued that unaccept- 
ably large biases, 6ioo — 5, would be required in order 
to make scaling defect theories with ^l„i — 1 compatible 
with current observations. We named this the 6ioo prob- 
lem and suggested that its resolution is the most serious 
challenge facing defect models for structure formation. 

These recent calculations (see also ref. ||]) represent 
the first calculations using a linear Einstein-Boltzmann 
solvei|j which includes all the relevant physical effects. 
The only input required is the unequal time-correlator 
(UETC) for the defect stress-energy, or an ensemble of 
source histories with the same UETC. In refs. we 
used a sophisticated model to represent the stress-energy 



^The present calculations and those of refs. used a mod- 
ified version CMBFAST which has been shown to give ac- 
curate (~ 1%) results for most inflationary models. 



of a network of cosmic strings, while numerical simula- 
tions were used in refs. ^,||. 

In the past, uncertainties in the results of defect cal- 
culations have meant that work has tended to concen- 
trate on the simplest possible models. Some attempts 
have been made to modify the semi-analytic approach 
of Albrecht and Stebbins to include additional pa- 
rameters, such as curvature (an open universe, 17^) or 
a cosmological constant (JIa) I^J^- This method yields 
only the shape of the matter power spectrum, although 
in ref. |^ a simple model was employed to estimate the 
dependence of the large angle fluctuation amplitude on 
fiA and f2fc, allowing rough statements about the rela- 
tive CMB/matter normalization to be made. The results 
of these calculations — reduced bias required to recon- 
cile theory with the data, and an improvement to the 
shape of the matter spectrum — serve as impetus for 
the current work, which considers a flat universe with 
a non- vanishing cosmological constant. Our calculations 
use the same modified version of CMBFAST and source 
two-point functions as in refs. and represent a self 
consistent calculation of both the CDM and CMB power 
spectra over all observable scales. We should note that 
similar results for the matter power spectra have been 
found by Avelino et al | |To[ |. 

The introduction of a cosmological constant will be 
seen to improve the required value of the bias in two 
ways. The first is in terms of the modified evolution of 
the background spacetime, which is illustrated in Fig. 
This shows the results for our standard scaling string 
source (described in detail in Q) with various choices of 
the parameter IIa, the proportion of the critical density 
from the cosmological constant with the Hubble constant 
Hq = 100/ikmsec~i Mpc"^ given by /i = 0.5. We find 
that as f2A is increased, the amount of power on large 
scales (around 100/i~^Mpc) is increased, with a corre- 
sponding decrease on small scales, bringing the shape of 
the matter spectrum in line with the galaxy data and 
reducing the magnitude of the feioo problem. This mod- 
ification can be understood via the standard shape pa- 
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FIG. 1. The (COBE normalized) matter power from the 
scahng cosmic string model for various values Qa- =0.0 
(short-long dashed line), 0.3 (dotted dashed line), 0.5 (short 
dashed line), 0.7 (long dashed line) and 0.9 (solid line). In- 
cluded also are the observed data points and the Standard 
Cold Dark Matter model (dot-dashed curve). 

rametcr F « Qmh and the equivalent shift in the time of 
radiation-matter equality. We find that Q\ = 0.7 seems 
to reproduce the shape of the observed spectrum partic- 
ularly well, although assuming perfect scaling this still 
requires a bias 6 w 4, over all scales. Note that we plot 
P{k)n'^^ rather than just P{k), so as to undo the redshift 
space distortion calculated in ref. JTl]] . 

In our previous work [^|^] we have made an extensive 
study of the effects of deviations from scaling. In the 
case of r^A = 0.0, we found that only very radical devia- 
tions from scaling can improve the amount of power on 
100/i~"'^Mpc scales, at the expense of creating an over pro- 
duction of power on smaller scales. However, given the 
improvement in the shape and amplitude of the matter 
power spectrum from the inclusion of a cosmological con- 
stant, one might think that mild deviations from scaling 
might improve things further. 

Rather than using the arbitrary scaling deviations to 
deduce a transition which yields the best case scenario, 
for the bulk of this work we have used the velocity de- 
pendent one-scale model for string evolution which 
was also used as the basis for the string evolution in ref. 
1^ . This allows us to calculate the three quantities which 
are variables in our model, namely the string coherence 
length relative to the horizon ^(t), the rms string velocity 
v{t) and the string density Ps(r) = /i/^^, as functions of 
conformal time, where /i is the mass per unit length of the 
string. The evolution equations for this model become 
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FIG. 2. The results of the velocity dependent one-scale 
model for the same values of f^A as in Fig. |l|. On left the 
density of strings relative to the background against a/oeq 
(bottom) and the rms velocity of the strin gs (top). On the 
right the same quantities but this time against r. 



is the chopping efficiency and fc(r) is the curvature of the 
strings. In order to interpolate between the measured 
values in the radiation and matter dominated eras, we 



use 
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where t is the length scale on the network, H — a' /a, c{t) 



where c,. = 0.23, k, = 0.17 and = 0.18, fcm = 0.49, 
the values in the radiation and matter eras respectively, 
and g = 0.19 is a constant which determines the speed of 
the transition. It has been suggested |l^ that this model 
reproduces the behaviour of a string network through the 
radiation-matter transition. 

The solution of these equations is illustrated in Fig. |^ 
for various values of fl^ and h = 0.5. We see that the 
density of strings shifts by approximately a factor of two 
in the interval a/acq=l to 100 independent of the value 
of S^A and in the models with Q\ ^ there is a fur- 
ther shift in the string density, during the A dominated 
phase; the effect being most obvious for large But 
when plotted against r, we see that r^q, and hence the 
transition era, is shifted as 17"^ in conformal time. For 
flm = 1, tbe velocity interpolates between the measured 
values in the radiation and matter eras, with a transient 
increase during the transition, and for Q\ =/= 0, the veloc- 
ity of the strings is exponentially suppressed during the 
A dominated era. Physically, the cosmological constant 
is causing the universe to expand exponentially, reducing 
the velocity of the strings, preventing reconnection and 
hence the strings are no longer scaling. 

Here, we comment briefly on the implementation of 
these deviations from scaling into our string model. The 
velocity and length of each string segment at each time 
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FIG. 3. The (COBE) normalized matter power spectrum 
for the same values of f^A as in Fig. but this time using 
the velocity dependent one scale model to calculate ^ and v. 
Once again we include the observed data point s and Standard 
CDM. 
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FIG. 4. The CMB power angular power spectrum for the 
same values of Qa as in Fig. hi. Also included are the ob- 
servations data points and the standard CDM curve, 
(dot-dashed) curve for comparison. 



are simply chosen to be v{t) and ^(r), but for techni- 
cal reasons we now pick a fixed value of v rather than 
choosing it from a random distribution, which has been 
shown to make very httle difference to the resulting power 
spectrum. The number density of strings is taken to 
be (?(T)r)^'^, that is, there is one string per correlation 
length cubed. 

Figs. H and ^ show the matter power spectra and CMB 
anisotropics for the same values of fl^ as Fig. |^, with ^, 
V and number density of strings given by this model. 
Coupled with the effects of the background, the devia- 
tions from scaling improve &ioo significantly for f^A 7^ 0, 
with the best fit for Ha = 0.7 (and hence, for a similar 
value of flm as favoured in refs. although the inclu- 

sion of the time dependence on ^ has slightly affected the 
shape of the power spectrum. We estimate that a bias of 
&100 = 2.6 would give a good fit to the observations on 
large scales at the expense of an overproduction of power 
on smaller scales. A robust feature of the CMB spectra 
appears to be a peak at around I = 400 — 600, whose 
height is increased as ft a increases. Finally, we should 
note that the value of string density per unit length G/i 
required to normalize to COBE is around I x 10~^, al- 
most independent of Q\, if we assume that the effective 
mass per unit length is /2 w I.7/i. This is well below the 
upper bound imposed by the absence of timing residuals 
in measurements of milli-second pulsars |jl^ . 

We have shown, therefore, that the inclusion of a cos- 
mological constant can improve the amount of matter 
power of larger scales, and that an interesting by prod- 
uct of this is a peak in the CMB angular power spec- 
trum. We now illustrate that once the lack of large scale 
power has been remedied, it is relatively simple improve 



the fit to the observations on smaller scales. Ignoring 
for the moment the constraints on n^h^ from Big Bang 
Nucleosynthesis (BBN), Fig. ^ contains three curves cre- 
ated using the velocity dependent one-scale model and 
riA = 0.7 with (I) Hb = 0.05 and h = 0.7, (2) Qb ^ 0.15 
and h = 0.7, (3) Hb = 0.125 and h = 0.5. Model (I) has 
the distinction that it fits the amplitude of the small scale 
matter data {bs ^ 1.0) though the shape is clearly not 
right, with a bias of 2.4 on 100/i-^Mpc scales. Models (2) 
and (3) rectify the small scale slope of the spectrum by 
increasing the baryon fraction, which damps out small 
scale power. A similar effect could be obtained by the 
introduction of hot dark matter. 

As one final point we return to the approach adopted 
in refs. |^,^ of finding the transition which fits the data 
best, since the exact nature of the deviation from scaling 
at the onset of cosmological constant domination has not 
yet been simulated. We implement a transition using the 
model described in ref. by varying /i with tt — 10000, 
Lt = 0.1 and x = 4, such that the time and length of 
the transition closely mirror those obtained in the one- 
scale velocity model and the amplitude chosen to give the 
best fit to the data. The results are shown in Fig. || first 
without bias and, to emphasize the goodness of the fit to 
the shape of the matter data, with a scale invariant bias 
of 6 = 2.0. This model has been shown by Gawiser and 
Silk ||l^ to give a good fit to the entire dataset, including 
CMB and matter power spectra, peculiar velocities and 
cluster abundances. 

Our results demonstrate that some of the problems 
of standard defect scenarios in underproducing power in 
the large scale matter distribution can be remedied by 
the introduction of a cosmological constant. Within this 
scenario, there are two distinct effects, a change in the 
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FIG. 5. The matter power spectrum for various models 
with Qa = 0.7. Velocity one scale model: (1) h = 0.7, 
Qb = 0.05 - long-short dash line; (2) h = 0.7, Qt = 0.15 
- short dash line; (3) h = 0.5, = 0.125 - dot-short dash 
line. Transition model with x = 4: unbiased - dotted curve; 
biased - solid curve. 



evolution of the background perturbations, and more im- 
portantly a modification to the scaling behaviour of the 
strings, which both act to allow more matter power on 
scales of 100/i~^Mpc. Using the velocity dependent one 
scale model we find that biases b ^ 2 are required: On 
first examination, this value still seems unreasonably high 
p6t , but the discrepancy with the data is small enough 
that these models deserve further consideration. The im- 
provement to the bias does not contradict the findings of 
refs. since it has come about by significant devia- 

tion in the scaling behaviour of the strings, due to the 
presence of the cosmological constant. Our results for the 
matter power spectrum and bias are in good agreement 
with ref. |^ , lending weight to the simpler techniques em- 
ployed in that work. It is clear that similar results will 
also be obtained in an open universe |g|jlo[|, and we are 
currently investigating this case. 

Finally, we should comment on our choice of cosmolog- 
ical parameters since this is a subject which is currently 
evolving rapidly. It appears that the best fit to the power 
spectrum is attained by using Qj^ = 0.7 which is very 
much in line with the current measurements of cosmolog- 
ical deceleration from supernova type la [ pT| and other 
independent measurements of the matter density of the 
universe. We should also note that some of the values 
of fib used in Fig. ^ are slightly larger than usually used 
on the basis of BBN. However, these higher values are 
not excluded ||l^ and we only use them here to illustrate 
that models which have a reasonable 6ioo can be made 
to fit the observed power spectrum on smaller scales — 
for example, the addition of hot dark matter can have a 
similar effect. We conclude that the interaction between 



this model and observations in the coming months will be 
interesting, and could in the end produce strong support 
for our model. 
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